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We study a retarded potential solution of a massless scalar field in curved space-time. In a special 
ansatz for a particle at rest whose magnitude of the (scalar) charge is changing with time, we found 
an exact analytic solution. The solution indicates that the phase velocity of the retarded potential of 
a non-moving scalar charge is position dependent, and may easily be greater than the speed of light 
at a given point. In the case of the Schwarzschild space-time, at the horizon, the phase velocity 
becomes infinitely faster than the coordinate speed of light at that point. Superluminal phase 
velocity is relatively common phenomenon, with the the phase velocity of the massive Klein-Gordon 
field as the best known example. We discuss why it is possible to have modes with superluminal 
phase velocity even for a massless field. 



PACS numbers: 



I. INTRODUCTION 



The speed at which fields (e.g. scalar, vector and grav- 
itational) propagate is a very subtle question. According 
to the special relativity, energy (and mass) can not travel 
with speed greater than the speed of light, but there is 
nothing restricting the speed of auxiliary fields like poten- 
tials. Potentials describe interactions, and interactions 
are mediated by virtual particles, which are off-shell and 
do not have any a priori preferred speed. 

For example, one is tempted to assume that the field of 
the non-moving source is frozen and does not propagate, 
until the source/charge moves and the field re- arranges 
its distribution. This would effectively mean that a static 
field is infinitely rigid and propagates with infinite speed. 
This may make sense classically, however, quantum me- 
chanically interactions arc fluctuations in space induced 
by virtual particles. Therefore, the situation is dynami- 
cal. To find the speed at which some interaction propa- 
gates, one has to calculate explicitly the effects of retar- 
dation in the Green's function of the field that mediates 
that interaction. 

It is well known that the QED vacuum structure can 
affect the propagation of light even in fiat space. The so- 
called Scharnhorst effect is a phenomenon in which light 
signals travel faster in between the two closely spaced 
conducting plates, than outside of the plates The 
reason is the Casimir effect, i.e. the vacuum polarization 
effect is suppressed in between the plates, so the photon 
loses less time propagating in between the plates than 
outside. This gives a hint that a massless particles do 
not always propagate at the speed of light in vacuum. 

There is even more counter-intuitive example in curved 
space-time. Namely, Drummond and Hathrell demon- 
strated in Q that vacuum polarization is sensitive to the 
curvature of spacetime. For example, for a photon propa- 



gating in a curved space, vacuum polarization can induce 
a modification of the wave equation in such a way so that 
in some cases photons travel at speeds greater than unity. 
The effect seems to be dispersive, and the phase velocity 
approaches the speed of light at high frequencies. Since 
the high-frequency limit of the phase velocity determines 
causality, it seems like causality is preserved in case. Ex- 
tensive discussion of this effect can be found in 0-Q]- 

These examples imply that propagation of quantum 
fields in curved space-time is a very subtle question, with 
many potential surprises. 

The simplest and perhaps the most instructive case to 
study will be the case of the scalar field. The reason 
is that the scalar field potential is not gauge dependent. 
The only freedom we have is to add an extra constant, i.e. 
ip — > ?A+const, which in turn has no dynamical effect and 
can be fixed by setting ip = at infinity. The best way to 
find out the propagation velocity is perhaps to study the 
Green's function of a field. Once the Green's function is 
found, we can analyze the retarded potential for a given 
field and infer the speed at which the signal propagates 
from a point to a point. However, the difficulty of finding 
the general Green's function for a field in a curved space- 
time makes this approach very difficult. Fortunately, the 
full space Green's function is not absolutely necessary 
to study the propagation phenomena. A case in which 
an observer observes a modulated source will be suffi- 
cient to study. Therefore, we will consider a massless 
scalar field potential for a stationary (non-moving) but 
time-dependent source. Our result shows that the phase 
velocity of the retarded potentials is position dependent, 
and may easily be faster than the speed of light. In the 
case of the Schwarzschild space, this phase velocity at 
the horizon can even be infinitely greater than the speed 
of light at the horizon. Though our solution in the lin- 
ear ansatz is analytic, our analysis of the general form 
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of the source is numerical. We therefore do not have an 
analytic form for a complete Green's function for an ar- 
bitrary source. However, the fact that the phase velocity 
of the scalar field varies locally is important. Among the 
other things it implies that gravity must affect the path 
of the massless scalar field, which for example should lead 
to the gravitational lensing effect for the massless scalar 
field. 

We would like to emphasis at the very beginning that 
throughout the paper we will use the term "signal" in 
a loose sense. We will call any change in the field a 
"signal" . While a non-moving source does not emit any 
real particles, the phase of the field will change, and the 
speed of that change (phase velocity) we will call the 
speed of the signal. While this is not a real signal or 
information in a strict sense (i.e. transmitted by the 
group velocity) it will have some important consequences. 



II. RETARDED GREEN'S FUNCTION FOR A 
MASSLESS SCALAR FIELD IN A CURVED 
SPACE-TIME 

As a referent point, we first show the retarded Green's 
function for a massless scalar field in the fiat space-time. 
The speed at which the signal propagates through space 
can be read off the retarded solution. Consider a point 
particle in Minkowski space carrying a massless-scalar- 
field charge at the origin. Let the magnitude of its charge 
increase (or decrease) in time as g(t). The equation of 
motion is 

- &J> - dfa - fify = 4nS(x)g(t) (1) 
The solution for the function ip is 







g(t-M) 

\r\ ' 
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where r = (x,y,z). The scalar field potential falls off 
with distance in flat space as 1/r. From the numerator, 
we sec that the signal travels from a point to a point with 
the speed of light c, i.e., if we increase the magnitude of 
the charge at the origin, the potential at a point r will 
be affected after time t = \ f\. So, it takes some time for 
a signal to propagate even if the source is static. This 
is best understood in terms of virtual particles. A static 
source emits the sea of virtual particles which modify the 
space around it. Result in Eq. © implies that, in flat 
space virtual (just like real massless) particles propagate 
with the speed of light, at least as long as they are in 
vacuum. 



III. GREEN'S FUNCTION FOR A MASSLESS 
SCALAR FIELD IN A CURVED SPACE-TIME 



There are very few examples of exact Green functions 
in curved space-times [8l-fl4j . The reason is that it is no- 
toriously difficult to find an exact solution without any 



approximations flM [l6j. However, we will demonstrate 
that the particular case with a great degree of symmetry, 
i.e. a charge located at the center of a spherical symmet- 
ric curved space, is directly solvable. We will then use 
the explicit solution to discuss the speed at which scalar 
field potentials propagate in such space-times. 

We fix again a point scalar charge at the origin of a 
spherically symmetric space. Fixing the charge at the ori- 
gin rather than at an arbitrary point in space will provide 
the required symmetry and greatly facilitate the problem. 
We let its magnitude change as git). The geometry of the 
space-time can be written as 



dr 2 



g tt dt 2 + g rr dr 2 - r 2 (dQ 2 + : 



(3) 



The equation that we will try to solve is 



D t D t ip + D r D r i\> + D e D e ip 



D*D^ = 6(r)^ (4) 



where h = —gugrrf 4 . 

We note here that our definition of the scalar charge 
slightly differs from the definition in Eq. 2.2 in 13j. The 
difference is the time component of the four- velocity u*, 
which in our static case (charge is not moving) is just a 
constant and can be absorbed in g(t). Further, strictly 
speaking, we are dealing with geometries without hori- 
zons, so we will not discuss the no-hair theorems. 

We will first try to find the time-independent solution, 
which will correspond to a static scalar charge of constant 
magnitude. In that case, Eq. Q reduces to 



l j=d r (g rr Vhd r ^) =S(r) 



(5) 



Without loss of generality, we can set g(t) = 1. The 
solution can be found by applying volume integration 
over the element y/hdtd 3 x. The static solution ip a is 



tps{r) 




-9tt(R) 1 

gu(R) R 2 



dR 



(6) 



It is easy to verify that this is a solution by substituting 
Eq. ([5]) back into Eq. Q . Since this solution depends on 
both g t t and g rr , it is different from the result from the 
flat space, but it will reduce to the flat space solution at 
large radius r. 

Wc will now try to find the time-dependent solution 
which corresponds to g(t). Since we a priory expect mass- 
less particle to propagate with the speed of light, we may 
expect this solution to have the following form 




Mr), (7) 



where the expression in parentheses on the right-hand 
side is the argument of the function g. This form is a 
straightforward curved space generalization of the flat 
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space solution given by Eq. @. The term w is just 



the coordinate speed of light in the radial direction (ob- 
tained from dr = 0). Though this appears to be a reason- 
able guess, this form is a solution only if g rr = —gtt- How- 
ever, this requirement brings us back to the flat space. 
We will therefore generalize the form of the solution al- 
lowing for the possibility that the propagation speed is 
not the speed of light. We now try a more general form 



i>(t, r )=g[t- 



1 



v(R) 



dR Mr) 



(8) 



where v(r) is the coordinate speed at which the signal of 
the retarded potential travels (not necessarily the speed 
of light). This v(r) must asymptotically go into the speed 
of light at large r where the space-time becomes flat. We 
do not expect this form to always generate a solution. 
However, if g is a linear function of its argument, one 
can find the suitable solution. Therefore we consider the 
linear form of g 



i 



v(R) 



dR 



A t 



1 



v(R) 



dR 



B (9) 



where A and B are two constants. If wc plug this ansatz 
into Eq. (gj, we find a condition under which the solution 
is valid 



~9tt 



2/2 i i A 

r ip s =r ip s c t 



2„/,2, 



(10) 



Here, q = J is the coordinate speed of light. Since 

v is also a coordinate velocity of propagation, it will be 
different for different observers and it will change from 
point to point. But, it is clear that the retarded signal 
propagates at a speed that is different from the speed of 
light for a given observer. 

To make this more clear, we consider the Schwarzschild 
space-time, i.e. 



9tt 



-g rr 



1 



2m 

r 



(11) 



We plug this condition into Eq. ([Id]), and the static solu- 
tion becomes 



1 



1 



1 



2m p2 
R Jl 



dR 



-m(l - ±f 
2m 



1 



+ - (12) 



We see that the static scalar field potential does not fall 
off as 1/r which was the case in flat space. However, in 
the limit of r 3> 2m, we recover the usual 1/r behavior. 

From Eq. (|TD)) . the coordinate propagation speed of the 
retarded potential is 



r 2 ip 2 ci > ct 



(13) 



which is not equal to the coordinate speed of light q, 
and in fact is always greater than c/ (for this particular 



example of the Schwarzschild space-time). Our calcu- 
lations will be strictly valid as long as our space-time is 
not strictly a black hole, but the conclusions will be valid 
even when we are only slightly outside the horizon. In 
the extreme limit, exactly at the horizon, r = 2m, the 
propagation speed, v, becomes infinitely faster than the 
speed of light at that point. 

In the context of the Schwarzschild black hole, the co- 
ordinate speed of light, q, vanishes at the horizon, and 
any signal sent from the horizon gets infinitely redshifted. 
Thus, in the standard description it remains unclear how 
information about the black hole charge which is presum- 
ably imprinted at the horizon can be communicated to 
the region around the black hole. [If the scalar charge 
is conserved, then formation of the black hole can not 
violate this conservation [l7[ •] Let's check what happens 
when the propagation speed of the retarded potential is 
taken into account. We can calculate the time, At, for a 
signal to propagate from the horizon to some finite dis- 
tance R 



At 



1 

—i 

2m v 
R 



-dr 



1 



dr 



nl-f [ r]n (l-^)]' 
-2m 



ln(l - 2m/ R) 



(14) 



This time is finite and therefore the potential has no 
problem to propagate from the horizon outside. Thus, 
a charged particle (at least with the scalar charge) can 
keep communicating its potential to the region outside 
the black hole. We mention again that our results are 
not strictly applicable to the black hole case, but we can 
always consider a shell whose radius is just slightly out- 
side its own Schwarzschild radius and preserve the qual- 
itative conclusions drawn here. In the next section we 
will reveal that the propagating velocity, v, is the phase 
velocity. Thus, this velocity refers to the change of phase, 
and it is not a group velocity. The casual light cone for 
real particles remains the same, the Green's function does 
not have support outside the light cone, and causality is 
preserved. 



IV. GENERAL CASE 

The discussion of the time dependent source so far 
was based on a solution found in the particular ansatz 
of Eq. ([9]) . We will now try to analyze the general form 
of g(t) (not only the linear ansatz that we used). We will 
first decompose the source g(t) into different frequency 
modes as 



9(t) = / fiKk-O exp(io;i)c?a; 



(15) 



The wave number of each frequency mode is oj/v^ir), 
where w w is the phase velocity of that mode. Then, the 
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scalar potential can be written as 



= / ff(w)exp 



iuj(t ■ 



-dr) 



f u (r)dij (16) 



where f u is the amplitude of the mode labeled by the 
frequency uj. g{ui) and g can be found by matching the 
boundary condition at r = 0. By plugging the above 
equation into equation we find that f u and v u must 
satisfy 



implies that becomes the speed of light q. It is this 
feature that ensures causality. 

For the to ^ modes, we will again use the spheri- 
cally symmetric Schwarzschild geometry. The boundary 
conditions are 



r ->• oo, P w = 1 
r — > oo, d r P u = 

^f = (i-iA) 2 



(23) 
(24) 
(25) 



9tt 



—7=O r ( Orfi 



Vu>g r : 



o fu + 
Vtgrr Vh V g r 







(17) 



except at r = 0. Since both imaginary and real parts 
must vanish independently, the above equation can be 
rewritten as two equations 



Jul 

9tt 



vigrr vh v g rr 



drfu 1 „ / Vh \ 

+ ~7r dr 

V u grr Vh \grrVu ' 



= 



(19) 



Eq. (jl9j) can be easily solved by integrating with re- 
spect to r. 



2 Vh 

1 



constant 



(20) 



In the limit of r — > oo, the space becomes flat, which 
implies v u — > 1 and f u — > 1/r, as it should. 
Eq. (f2T))) can be rewritten as 



P: 



-gu 



9r 



(21) 



where P u = f u r. If we substitute this relation into 
Eq. (TT5)) , and replace f u r with P u we get 



d^P u + -8 r ln 



-gu 

grr 



d r Pu - Tjdr hi 



~9tt 

grr 



Pu 

r 



(22) 



The zero mode, u> = 0, solution to this equation 
is exactly f u = t/j s , where tp s is the time-independent 
g =const solution given in Eq. ([6]). Moreover, the phase 
velocity v u in this case is the same as the propagation ve- 
locity given by Eq. (fT0|) in the ansatz solution we found. 
This then reveals the meaning of the parameter v in Sec. 

inn 

In the high frequency limit, u> — ¥ oo, the last term 
in Eq. (f2"2"]l dominates. In order to satisfy the equation 
it has to vanish, thus requiring P u = 1. Eq. (|21l) then 
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FIG. 1: This figure shows Pu>( r ) for three different values of w, 
i.e. lu — 1Q _1 , uj — 1CF 2 and uj = 1CP 3 . We see that P w grows 
as it is approaching the origin, higher ui modes increase slower 
than lower ui modes, and > 1 everywhere. This behavior 
implies that the phase velocity u u is always superluminal (for 
this case of the Schwarzschild geometry), and that lower uj 
modes propagate faster than the higher uj modes. 

In Fig. [TJ we show P w as a function of r for several 
values of the frequency oj. We can see that P u grows 
as it is approaching the origin. It is also apparent that 
higher uj modes increase slower than lower ui modes. We 
do not plot P u near the horizon, because the singularity 
will cause numerical instabilities. Since P u > 1 for all 
values of r, the phase velocity defined by Eq. (|2"Tj) 



-gu 



9r 



> Cl 



(26) 



is greater than the speed of light everywhere. 



V. COMMON FEATURES WITH OTHER 
EXAMPLES WITH SUPERLUMINAL PHASE 
VELOCITY 

In this section we discuss some other known examples 
where the phase velocity is superluminal, which may have 
have something in common with our results. 
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Perhaps the best known example is that of a massive 
Klein-Gordon field in flat space-time. The dispersion re- 
lation is simply lo 2 = k 2 + to 2 . The phase velocity 



uj/k 



(27) 



is always greater than unity as long as to ^ 0. Moreover, 
superluminality is most pronounced for low frequencies 
ui, while for large frequencies we have u> « k, i.e. v p rs 1. 
However, the group velocity 



v g — dio/dk 



(28) 



is always less than unity. Comparing this result with the 
results we obtained in curved space, we might conclude 
that the curvature of space induces an effective mass to 
the massless scalar field, making it formally equivalent to 
the massive Klein- Gordon field with supcrluminal phase 
velocity. 

The other, less known example is a massless scalar field 
in a (5 + l)-dimensional flat space-time. The wave equa- 
tion is 

d^-dl^-dl^-dl^-dl^-dl^ = s(x)S(t) (29) 

The general solution can be found in most mathematical 
physics textbooks (e.g. [liij ) or papers [l9| . The Green's 
function for this case is 



where, r = \J x\ + x\ + x\ + - 1 ' 2 
the following concrete source 



3 -r x\ + Xg. If one considers 



f(t,x) =sm(ut)5(x) (31) 
then the wave function can be easily found as 

ip = J f(t', x')G 5+1 (t-t',x- x>)dt'dx> (32) 
1 (bj cos(o;(t — r)) sin(cj(t — r)) ' 



8tt 2 



(33) 



Since this form includes two trigonometric functions, it is 
hard to see how the phase changes. We will then combine 
the two terms into a single trigonometric function. 



ip = S(u>,r)sm(u(t-R)) 



S{U},r) = -—J ; 

R = r — 4>{r)/bj 



( r = sm 



VI + r 2 uj' 



(34) 

(35) 
(36) 
(37) 



This form is similar to the form we studied in the last 
section. The phase velocity in this case is 



1 1 
1 



p d r R~^ [ r 2 co 2 



(38) 



We can see that the phase velocity is infinite at the 
origin (for fixed ui) but equal to the speed of light at r — > 
oo. In this case the solution is created by two waves with 
different phases. Since their amplitudes decay in different 
ways, their combination makes the total phase velocity 
change with location, and in fact makes it infinitely faster 
than the speed of light at some locations. These are 
the features which are shared with our solution for the 
massless scalar field in a curved space. Moreover, for a 
fixed finite r, superluminality is again most pronounced 
for small oj, while for large frequencies we have v p ~ 1. 



VI. CONCLUSIONS 

In this paper we analyzed the question of the speed at 
which potentials propagate in curved space-time. While 
finding an answer is easy in flat space, it becomes highly 
non-trivial in curved space-time. The difficulties range 
from finding an exact solution for the Green's function 
to choosing the right definition of the propagation speed. 
To avoid gauge and other ambiguities we considered the 
massless scalar potential. We located the scalar charge 
whose magnitude was changing in time at the origin in 
a spherically symmetric space-time, and found the so- 
lution for different frequency modes for this configura- 
tion. A non-moving particle does not emit real scalar 
field quanta, but what is changing in the system is the 
phase of the field. We found that the phase velocity is 
not constant but changes from point to point. Moreover, 
in the specific case of the Schwarzschild geometry, it is 
always greater than the coordinate speed of light at any 
given point. In an extreme limit, exactly at the horizon, 
the phase velocity becomes infinity faster than the speed 
of light at that point (which is actually vanishing). In 
fact, this feature is required if a black hole is going to 
communicate "information" about its potential which is 
presumably located at the horizon to the outer world. 

It is important to note that the phase velocities, v u , 
for different frequency modes (labeled by the frequency 
oj) are different for each mode, and in general they are 
different from the local speed of light. Also, the am- 
plitudes of different frequency mode's (f u in the text), 
have different r dependence. These two facts make the 
curved space case quite different from the flat space, and 
explain why it was impossible to find a uniform prop- 
agation mode like in flat space (see Eq. ©. Since the 
propagation is dispersive, the high-frequency limit of the 
phase velocity will determine causality. Since the phase 
velocity approaches the speed of light at high frequencies, 
causality is preserved in our case. 

The cases of the electromagnetic and gravitational po- 
tentials are more complicated because of the non-zero 
spin. However, they are also massless fields and will per- 
haps have some similar properties. In particular, we ex- 
pect the retarded electromagnetic and gravitational po- 
tential from a non-moving source to propagate at a speed 
different from the speed of light. In other words, the aver- 
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age velocities of virtual photons and gravitons should not 
be the same as for real photons and gravitons in curved 
space-times. 

A related question can be asked in the context of gravi- 
tational lensing. If the retarded gravitational potential of 
a static source travels with a finite speed (not necessarily 
the speed of light), it must experience the effect of the 
gravitational lensing, just as the light does. This would 
imply that the gravitational lensing effect on gravitons 
should be able to amplify or reduce the strength of grav- 
ity from a given static source [2(| ■ In [2l[ , several exam- 
ples were constructed to emphasis that the gravitational 
lensing could affect real gravitons, but could not lens any 
static gravitational field potential (though the static po- 
tential could be affected to some extent). However, the 
sources used in these examples were infinite planes, and 
not point sources, so the conclusions are perhaps not gen- 
eral. 

If our conclusions for the scalar field hold for gravi- 
tons as well, then the static gravitational potential could 
propagate at any finite speed (except in the extreme case 
of the black hole horizon where it should be infinite) de- 
pending on the curved background. Since this speed is 
finite and position dependent, the effect of gravitational 
lensing of gravity should exist, though the magnitude of 
the effect should be different from the gravitational lens- 
ing of the light because of the different speed of propa- 
gation. It is interesting that one of the possible expla- 
nations of the Pioneer anomaly [22| is the focusing of 
gravity at around 25 AU [2(J, exactly where the Pioneer 
anomaly arises. This could be a hint that gravity is bent 
nearby our Sun [23| . of course if the real explanation is 
not something more conventional, like the thermal radi- 
ation pressure [H, [25| . 

At the end we would like to compare our findings with 
the existing similar results in the literature, e.g. In 
[l|, using the Casimir effect, the authors showed that vac- 
uum polarization effects may lead to supcrluminal prop- 
agation of photons in between the plates (since the vac- 
uum polarization effects arc suppressed there). While 
this is a flat space result, it is indicative that superlu- 



minality may arise in completely physical setups. In Q 
it was argued that the quantum corrections in curved 
space-time are able to introduce tidal gravitational forces 
on the photons which in general alter the characteristics 
of propagation, so that in some cases photons travel at 
speeds greater than unity. In that case it is actually the 
low-frequency limit of the phase velocity that is supcrlu- 
minal. This indicates that propagation of quantum fields 
in curved space-time is a very non-trivial problem, and 
surprising results may be derived. It should be noted that 
super luminality does not always lead to paradoxes, since 
in both of the above mentioned cases it is impossible to 
send signals backward in time. While work presented in 
[H is perturbative, our analysis is exact since it based 
on the exact solution of the Green's function in curved 
space-time. It is interesting that our analysis also indi- 
cate that lower frequency modes propagate faster than 
high frequency modes, in good agreement with Q. 

Finally, we emphasis again that we found only the 
phase velocity to be supcrluminal. If the group veloc- 
ity is not superluminal, then the Green's function does 
not have support outside the light cone, and causality is 
preserved. Strictly speaking, a second order linear wave 
equation can not have a "wavefront" propagating faster 
than the speed of light. However, this statement does not 
affect velocity of an individual frequency component of 
the phase. It is only when one takes into account all the 
frequencies (where the higher frequencies give the domi- 
nant contribution) that he has to obey that statement. 
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